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Stability of the Period-Doubled Core of 
the 90° Partial in Silicon 

In a recent Letter [1], Lehto and Oberg (LO) investi- 
gated the effects of strain fields on the core structure of 
the 90° partial dislocation in silicon, especially the in- 
fluence of the choice of supercell periodic boundary con- 
ditions in theoretical simulations. Specifically, they ad- 
dressed the relative stability of the traditionally accepted 
single-period (SP) geometry vs. our recently proposed 
double-period (DP) structure [2]. Performing supercell 
calculations that employed the original Keating poten- 
tial for Si [3], they reached two main conclusions: (1) 
that a dipole-cell configuration minimizes the overall cell 
strain, thus being more adequate for small-cell simula- 
tions; (2) that the relative stability of the SP and DP ge- 
ometries depends on the choice of boundary conditions, 
with the SP and DP cores being favored for "dipole" and 
"quadrupole" configurations, respectively. The purpose 
of the present comment is not to dispute their first point, 
with which we agree. Rather, we wish to focus on the 
more important DP-versus-SP stability issue, because we 
believe their second conclusion to be incorrect. Below, we 
show that their results for the relative stability between 
the two structures are in disagreement with cell-size con- 
verged tight-binding total energy (TBTE) calculations, 
which suggest the DP core to be more stable, regardless 
of the choice of boundary condition. Moreover, we argue 
that this disagreement is due to their use of a Keating 
potential. 

Clearly, in the limit of sufficiently large supercells, all 
results should be independent of the the choice of bound- 
ary conditions. Here, supercell convergence is investi- 
gated by performing TBTE and Keating-potential cal- 
culations for three different supercell sizes, using both 
the dipole and quadrupole cells, and two different sets of 
parameters for the Keating-potential calculations [3,4]. 
In our larger cell, all dislocation separations are similar 
to those of the 2048- atom cells studied by LO. The first 
observation we can draw from the TBTE results in Ta- 
ble I is that, as pointed out by LO, for small cells the 
dipole boundary condition gives results which are closer 
to the converged value, while the quadrupole cell has a 
bias of ~10 meV/A in favor of the DP structure, which 
decreases by one order of magnitude as we approach cell- 
size convergence. However, Table I also shows that the 
Keating potential has a much stronger bias in favor of 
the SP structure. For the 192-atom cell, the results of 
LDA (from Ref. [2]), TETB, and the Keating potentials 
of Refs. [3] and [4] for Edp-Esp are -80, -70, -20, and 
— 15, respectively (in meV/A, with error bars of ~±5 in 
each case). Thus, the Keating potential has a systematic 
bias in favor of SP of ~50 meV/A relative to the more 
accurate methods. As it happens, going to cell-size con- 
vergence shifts all the Ejjp — Esp values by roughly 20 
meV/A. As a result, the converged values for the Keating 


potentials are close to zero, making the variations with 
choice of boundary conditions (which, in absolute num- 
bers, are similar to those seen in TBTE) take on an arti- 
ficial importance. In fact, the Eup — Esp values are so 
small that even the qualitative conclusions about which 
structure is favored can be seen to depend on the choice 
of Keating parameters. On the other hand, it seems clear 
that the more accurate TBTE results would be quite im- 
mune from displaying a sign change in Edp — Esp as 
a result of details of the choice of cell size or boundary 
conditions. 

Finally, LO also based their conclusions in part on a 
density-functional calculation for a finite cylindrical sam- 
ple. However, we suggest that the systematic errors in- 
volved in choosing surface boundary conditions for such 
a geometry may also be severe, and that such results 
should not be trusted in the absence of careful tests of 
convergence with respect to sample size. 

TABLE I. Energy of tlie DP relative to tlie SP core, in 
me v/A per dislocation, for the 90° partial in Si. Four dif- 
ferent approximations (LDA, TBTE and two different Keat- 
ing parameterizations) are used to compare quadrupole and 
dipole boundary conditions. Cell size refers to the DP case. 


Cell size (atoms) 

192 

576 

1920 

LDA 




quadrupole 

-79 



TETB 




quadrupole 

-74 

-63 

-55 

dipole 

-62 

-55 

-55 

Keating, Ref. [3] 




quadrupole 

-27 

-7 
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dipole 

-14 

2 

5 

Keating, Ref. [4] 




quadrupole 

-22 

-10 

-5 

dipole 

-13 

-5 

-3 
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